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ABSTRACT 

In this paper we consider stationary force-free magnetosphere of an aligned rotator 
when plasma in the open field line region rotates differentially due to presence of a 
zone with the accelerating electric field in the polar cap of pulsar. We study the impact 
of differential rotation on the current density distribution in the magnetosphere. Using 
split-monopole approximation we obtain analytical expressions for physical parameters 
of differentially rotating magnetosphere. We find the range of admitted current density 
distributions under the requirement that the potential drop in the polar cap is less 
than the vacuum potential drop. We show that the current density distribution could 
deviate significantly from the "classical" Michel distribution and could be made almost 
constant over the polar cap even when the potential drop in the accelerating zone is 
of the order of 10 per cents of the vacuum potential drop. We argue that differential 
rotation of the open magnetic field lines could play an important role in adjusting 
between the magnetosphere and the polar cap cascade zone and could affect the value 
of pulsar breaking index. 
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1 INTRODUCTION 

The physics of radio-pulsars is still not fully understood de- 
spite large effort of many theoreticians in this field. It is 
generally assumed that radiopulsar has MHD-like magne- 
tosphere which is very clos e to the force- free s tage - the 
model firstly introduced bv lGoldreich fc Julian! (|l969u . For 
many years solution of force-free MHD equations was the 
problem, even for the simplest case of an aligned pulsar. 
Now to solve the Grade-Shafranov equation describing struc- 
ture of a force-free magnetospher e of an aligned pulsar is 
not a problem anymore (see e.g. IContopoulos et aD 1 19991 ; 
iGoodwin et"aHl2004 lGruzinovll2005l ; iTimokhinl 120061 ). Sta- 
tionary magnetosphere configurations for an aligned rota- 
tor were obtained als o a the final stage in non-stat i onary 
numerical modelling (iKomissarovl 120061; iMcKinnevI 120061 ; 
iBucciantini et al.1 120061 ; ISpitkovskvl l2006h . Eve n the case 
of an inclined rotator was studied numerically l|Spitkovskv| 
2006). However, the pulsar magnetosphere is a very compli- 
cated physical system because most of the current carriers 
(electrons and positrons) are produced inside the system, 
in the polar cap cascades. Production of electron-positron 
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pairs is a process with a threshold, so it could operate only 
under specific conditions and, generally speaking, not any 
given current density could flow through the cascade zone. 

In magnetohydrodynamics (MHD) the current density 
distribution is not a free "parameter", it is obtained in 
course of solving of MHD equations. In case of pulsars ob- 
taining a solution of MHD equations does not solve the prob- 
lem, because it could happen that the polar cap cascade zone 
could not provide the required current density distribution 
and, hence, support the particular configurations of the mag- 
netosphere. In terms of MHD the polar cap cascade zone sets 
complicated boundary conditions at the foot points of the 
open magnetic field lines and any self-consistent solution of 
the problem must match them. The most "natural" configu- 
ration of the magnetosphere of an aligned rotator, when the 
last closed field line extends up to the light cylinder, requires 
current density distribution which could not be supported 
by stationa ry electromagnet ic cascades in the polar cap of 
pulsar (see iTimokhinl [20061 . hereafter Paper I). That con- 
figuration requires that in some parts of the polar cap the 
electric current flows against the preferable direction of the 
accelerating electric field. This seems to be not possible also 
for non-stationary cascades, although this problem requires 
more carefully investigation than it has been done before 
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i|Fawlevlll978l ; lAl'Ber et ai1ll975l ; ILevinson et alj|2003 ). So, 
the structure of the magnetosphere should be different from 
this simple picture. The magnetosphere of a pulsar would 
have a configuration with the current density distribution 
which can flow through the polar cap cascade zone without 
suppression of electron-positron pair creation. Whether such 
configuration exists is still an open question, i.e. a possibility 
that the real pulsar magnetosphere has large domains where 
MHD approximatio n is broken could not be completely ex- 
cluded too (see e.g. lArondll979l ; lMichellll99lh . 

As the pulsar magnetosphere and the polar cap cascade 
zone have too different characteristic timescales, it would 
be barely possible to proceed with modelling of the whole 
system at once. Therefore, these physical systems should be 
modelled separately and the whole set of solutions for each 
system should be found, in order to find compatible ones. 
Namely, we suggest the following approach to the construc- 
tion of the pulsar magnetosphere model: one should find out 
which currents could flow through the force-free pulsar mag- 
netosphere and compare them with the currents being able 
to flow through the polar cap cascade zone. In this work 
we deal with the first part of the the suggested "program" . 
Namely, we consider the range of possible current density 
distributions in force-free magnetosphere of an aligned rota- 
tor. 

Force-free magnetosphere of an aligned rotator is the 
simplest possible case of an MHD-like pulsar magnetosphere 
and needs to be investigated in the first place. This sys- 
tem has two physical degrees of freedoms i) the size of the 
closed field line zone, and ii) the distribution of the angu- 
lar velocity of open magnetic field lines. In each stationary 
configuration the current density distribution is fixed. Con- 
sidering different configurations by changing (i) and (ii) and 
keeping them in reasonable range the whole set of admitted 
current density distributions can be found. Differential ro- 
tation of the open field lines is caused by variation of the 
accelerating electric potential in the cascade zone across the 
polar cap. Theories of stationary polar cap cascades predict 
rather small potential drop and in this case only one degree 
of freedom is left - the size of the zone with closed magnetic 
field lines. This case was studied in details in Paper I, with 
the results that stationary polar cap cascades are incom- 
patible with stationary force-free magnetosphere. So, most 
probably the polar cap cascades operate in non-stationary 
regime. For non-stationary cascades the average potential 
drop in the accelerating zone could be larger than the drop 
maintained by stationary cascades. Hence, the open mag- 
netic field lines may rotate with significantly different angu- 
lar velocities even in magnetospheres of young pulsars. On 
the other hand, for old pulsars the potential drop in the 
cascade zone is large, and magnetospheres of such pulsars 
should rotate essentially differentially anyway. 

The case of differentially rotating pulsar magnetosphere 
was not investigated in details before Although some authors 
addressed the case when the open magnetics field lines rotate 
differently than the NS, but only the case of constant angular 
veloc ity was considered fe.g. lBeskin et a l. 1993; C ontopoulosl 
2005). The first attempt to construct a self-consistent model 
of pulsar magnetosphere with differe ntially rotating open 
field line zone was made in iTimokhinl l)2007l 'l. hereafter Pa- 
per II. In that paper we considered only the case when the 
angular velocity of the open field lines is less than the an- 



gular velocity of the NS. We have shown that the current 
density can be made almost constant over the polar cap, al- 
though on a cost of a large potential drop in the accelerating 
zone. The angular velocity distributions was chosen ad hoc 
and the analysis of the admitted range for current density 
distributions was not performed. 

In this paper we discuss properties of differentially ro- 
tating magnetosphere of an aligned rotator in general and 
elaborate the limits on the differential rotation. We study in 
detail the case when the current density in the polar cap is 
a linear function on the magnetic flux. It allows us to obtain 
main relations analytically. We find the range in which phys- 
ical parameters of the magnetosphere could vary, requiring 
that a) the potential drop in the polar cap is not greater 
that the vacuum potential drop and b) the current in the 
polar cap does not change its direction. 

The plan of the paper is as follows. In Section [2] we 
discuss basic properties of differentially rotating force-free 
magnetosphere of an aligned rotator and derive equations 
for angular velocity distribution, current density and the 
Goldreich- Julian charge density in the magnetosphere. In 
Section [3] we derive equations for the potential drop which 
supports configurations with linear current density distribu- 
tion in the polar cap of pulsar and give their general so- 
lutions. In Section [4] we analyse the physical properties of 
admitted magnetosphere configurations: the current density 
distribution, the maximum potential drop, the angular ve- 
locity of the open magnetic field lines, the Goldreich- Julian 
current density, the spindown rate and the total energy of 
the magnetosphere. At the end of that section we we con- 
sider as examples two sets of solutions: the one with constant 
current densities and the another one with the smallest po- 
tential drops. In Section[5]we summarise the results, discuss 
limitation of the used approximation and briefly describe 
possible modification of the obtained solutions which will 
arise in truly self-consistent model. In that section we also 
discuss the issue with the pulsar braking index. 



2 DIFFERENTIALLY ROTATING 

MAGNETOSPHERE: BASIC PROPERTIES 

2.1 Pulsar equation 

Here as in Papers 1,11 we consider magnetosphere of an 
aligned rotator that is at the coordinate origin and has dipo- 
lar magnetic field. We use normalisations similaiQto the ones 
in Paper I, but now we write all equations in the spherical 
coordinates (r,9,(f>). We normalise all distances to the light 
cylinder radius of the corotating magnetosphere Rlc = c/O, 
where f2 is the angular velocity of the neutron star (NS), c is 
the speed of light. For the considered axisymmetric case the 
magnetic field can be expressed through two dimensionless 
scalar functions and 5* as (cf. eq. (8) in Paper I) 

H V^xe + 5e^, 
Rl c rsintf ' 1 ' 

where e$ is the unit azimuthal, toroidal vector, fi = 
Bo-Rns/2 is the magnetic moment of the NS; Bo is the mag- 
netic field strength at the magnetic pole, Rns is the NS 

1 note that here in contrast to Paper I \E" is already dimensionless 
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radius. The scalar function $ is related to the magnetic flux 
as $ mag (W, Z) = 2n (fj,/Rhc) y(r,0). $ mag is the magnetic 
flux trough a circle of a radius w = r sin 9 with its cen- 
tre at the point on the rotation axis being at the distance 
Z = rcos9 from the NS. The lines of constant ^ coincides 
with magnetic field lines. The scalar function 5* is related 
to the total current J outflowing trough the same circle by 
J(m,Z) = 1/2 (fi/Rl c )cS(r,9). 

The electric field in the force-free magnetosphere is 
given by 



(2) 



where /3 is the ratio of the angular velocity of the magnetic 
field lines rotation Of to the angular velocity of the NS, 
= Of/O (cf. eq. (14) in Paper I). The difference of the 
angular velocity of a magnetic field line from the angular 
velocity of the NS is due to potential drop along that line in 
the polar cap acceleration zone. 

For these dimensionless functions the equation describ- 
ing the stationar y force-free ma g netosphere, the so-called 
pulsar equation dMichell Il973al ; IScharlemann fc Wagoner! 
1 19731 ; IOkamotolll974h . takes the form (cf. eq. 20 in Paper I) 



((3rsm9) 2 ] A* 



+ S 



cos 9 de* 
sin 6 r 



+ 



dS_ 



P^r <rsin9 V*) 2 =0. 



(3) 



This equation express the force balance across the magnetic 
field lines. At the light cylinder the coefficient by the Lapla- 
cian goes to zero and the pulsar equation reduces to 



dS _ 



(V*) 2 +2/? sine (<9 r * + /3cos<9de*) . (4) 

Each smooth solution must satisfy these two equations and 
the problem of a solving the pulsar equations transforms to 
an eigenfunction problem for the poloidal current function 
S (see e.g. Section 2.3 in Paper I). Equation Q could also 
be considered as an equation on the poloidal current. 

We adopt for the magnetosphere the configuration with 
the so-called Y-null point. Namely, we assume that the mag- 
netosphere is divided in two zones, the first one with closed 
magnetic field lines, which extend from the NS up to the neu- 
tral point being at distance xo from the NS, and the second 
one, where magnetic field lines are open and extend to infin- 
ity (see Fig. [TJ . In the closed magnetic field line zone plasma 
corotates with the NS, there is no poloidal current along 
field lines, and the magnetic field lines there are equipoten- 
tial. Apparently this zone can not extend beyond the light 
cylinder. In the the rest of the magnetosphere magnetic field 
lines are open due to poloidal current produced by outflow- 
ing charged particles. The return current, needed to keep 
the NS charge neutral flows in a thin region (current sheet) 
along the equatorial plane and then along the last open mag- 
netic field line. We assume that this picture is stationary on 
the time scale of the order of the period of the NS rotation. 
As it was outlined in Paper I, the polar cap cascades in pul- 
sars are most probably non-stationary. The characteristic 
time scale of the polar cap cascades ~ h/c ~ 3 • 10~ 5 sec (h 
is the length of the acceleration zone being of the order of 
-Rns) is much shorter that the pulsar period (for most pul- 
sars being S> 10 -3 sec). So, for the global magnetosphere 




Figure 1. Structure of the magnetosphere of an aligned rota- 
tor (schematic picture). Magnetic field lines are shown by solid 
lines. Outflowing current J along open magnetic field lines and 
returning current Jroturn in the current sheet, separating the open 
and closed magnetic field line zones, are indicated by arrows. The 
current sheet is along the last open magnetic field line, corre- 
sponding to the value of the flux function V&pc- Distances are 
measured in units of the light cylinder radius for the corotating 
magnetosphere i?LC> i- e - the point with x = 1 marks the posi- 
tion of the light cylinder in the corotating magnetosphere. The 
null point x$ could lie anywhere inside the interval [0, 1]. Possi- 
ble positions of the real light cylinder are shown by dotted lines. 
The line I corresponds to the case when l/f3(^/ pc ) < xq; II - to 
x < l/Z^pc) < 1; III - to 1 < l//3(>I'pc) (see the text further 
in the article). 



structure only time average of the physical parameters con- 
nected to the cascade zone are important. In the rest of the 
paper, when we discuss physical parameters set by the cas- 
cade zone we will always mean the average values of them, 
unless the opposite is explicitly stated. 

Differential rotation of the open magnetic field lines 
which is caused by presence of a zone with the accelerat- 
ing electric field in the polar cap of pulsar i) contributes to 
the force balance across magnetic field lines (the last term 
in eq.Q), ii) modifies the current density in the magneto- 
sphere (the first term in r.h.s. of eq. (H}), and iii) changes 
the position of the light cylinder, where condition @ must 
be satisfied. Note that for the case i) and ii) the deriva- 
tive dj3/d^ , i.e. the form of the distribution /3(^), plays an 
important role. So, for different angular velocity distribu- 
tions in the open magnetic field line zone there should exist 
different magnetosphere configurations that have in general 
distinct current density distributions. Let us now consider 
restrictions on the differential rotation rate /3(\t). 



2.2 Angular velocity of the open magnetic field 
lines 

Due to rotation of the NS a large potential difference arises 
between magnetic field line foot points at the surface of the 
NS. The potential difference between the pole and the mag- 
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netic field line corresponding to the value of the magnetic 
flux function ^ is 



(5) 



In perfectly force- free magnetosphere the magnetic field lines 
are equipotential. However, due to presence of the polar cap 
acceleration zone, where MHD conditions are broken a part 
of this potential difference appears as a potential drop be- 
tween the surface of the NS and the pair-formation front, 
above which the magnetic field line remains equipotential. 
This potential drop is the reason why the open magnetic field 
lines rotate differently from the NS. The normalised angular 
velocity of a magnetic field line is expre ssed tro u gh th e 
potential drop along the field line as (e.g. iBeskinl {2005), 
Paper I) 



(3 = 1 + 



Rtc 



(6) 



V is the total potential drop (in statvolts) along the mag- 
netic field line in the polar cap acceleration zone (cf. eq. (23) 
in Paper I). 

The polar cap of pulsar is limited by the magnetic field 
line corresponding to a value of the flux function "Jpc. The 
potential drop between the rotation axis and the boundary 
of the polar cap is 



AVn 



(7) 



This is the maximum available potential drop along an open 
magnetic field line. It could be achieved in vacuum, when 
there is no plasma in the polar cap. We will call AV pc the 
vacuum potential drop. Let us normalise the poloidal flux 
function ^ to its value at the last open magnetic field line 
^pc and introduce a new function ip = ^/^pc. Normalising 
potential drop along field lines to the vacuum potential drop 
and introducing the dimensionless function V = l> / A V pc we 
rewrite expression for the normalised angular velocity of of 
the open magnetic field line as 



(8) 



As the potential drop along any field line can not be greater 
than the vacuum drop and could not have different sign than 
the vacuum drop, variation of the electric potential across 
the polar cap can not exceed the vacuum potential drop. In 
terms of the dimensionless functions this condition has the 
form 



|V(^l)-V(^2)|<l, V^l,V2 £ [0,1]. 



(9) 



Inequality <(9j sets the limit on the electric potential in the 
polar cap of pulsar. 

2.3 Current density in the polar cap 

In order to obtain the current density distribution in the po- 
lar cap of pulsar the pulsar equation <(3j together with the 
condition at the light cylinder Q must be solved. There 
is an analytical solution of the pulsar equation only for 
split-monopole configuration of the poloidal magnetic field. 
Namely, when flux function $ has the form 



being a constant, equations © and ((4| have a smooth 
solution if the poloidal cur rent function S has the form (e.g. 
iBlandford fc Znaieldll977h 



= -/?(*) *(2-^-). 

Wm 



(11) 



Here corresponds to the value of the magnetic flux 
trough the upper hemisphere, i.e. it corresponds to the mag- 
netic field line lying in the equatorial plane. The poloidal 
current given by equation JTTJ i s very similar to current in 
the well known Michel solution (|Miche]||l973bh , but this ex- 
pression is valid for non-constant /3(>t) too. 

In this paper we will use for the poloidal current func- 
tion S expression (JTTJ . Doing so, we assume that in the 
neighbourhood of the light cylinder the geometry of the 
poloidal magnetic field is close to a split monopole. This 
is good approximation if the size of the closed magnetic 
field line zone is much smaller than the light cylinder size, 
%o <C l//3(t/0 5 ip < 1. For configurations where the size 
of the corotating zon^f] approaches the light cylinder the 
poloidal current S is different from the one given by eq. (|11[) . 
but we expect that this deviation should not exceed 10-20 
per cents. Indeed, in the numerical simulations described in 
Paper I, where the case of constant /3 = 1 was considered, 
the deviation of S from the Michel's poloidal current did 
not exceed 20 per cents and it got smaller for smaller size 
of the corotating zone (see Fig. 3 in Paper I). Similarly, in 
Paper II, where we considered the case of variable /3 < 1, 
the poloidal current deviated from the values given by the 
analytical formula |TT} by less than 20 per cents and the 
difference became smaller for smaller size of the corotating 
zone. So, we may hope that the same relation holds in the 
general case too. 

We intent to find the range of admitted current den- 
sity distributions in the force-free magnetosphere. Here we 
use the split-monopole approximation for the poloidal cur- 
rent l|ll[). hence, we can study only the effect of differential 
rotation on the current density distribution. The dependence 
of the current density on the size of the corotating zone in 
differentially rotating magnetosphere will be addressed in a 
subsequent paper, where we will refine our results by per- 
forming numerical simulations for different sizes of the coro- 
tating zone. 

So, in our approximation the last closed field line in 
dipole geometry corresponds to the field line lying in the 
equatorial plane in monopole geometry, i.e. = ^pc- In 
normalised variables expression for the poloidal current has 
the form 



S(*) = -* pc /3(V>M2-V). 



(12) 



The poloidal cur rent density in the magnetosphere is (see 
e.g. lBeskinll2005l) 

_ c (j, VS x e<f, _ fi-Bpd 1 dS . « 



4tt< c 



* = *m(1 - cos( 



(10) 



In the polar cap of pulsar the magnetic field is dipolar and, 
hence poloidal. The Goldreich- Julian charge density for the 



2 plasma in the closed field line zone corotates with the NS, so 
we will call the region with the closed magnetic field lines the 



corotating zone 
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corotating magnetosphere near the NS is 



o 

Pgj 



n b 

2-kc 



(14) 



Using expressions (|12p - (114p we get for the current density in 
the polar cap of pulsar 



ij GJ [2/3(1 -V0 +/3V(2 -V)] 



(15) 



The prime denotes differentiation with respect to ip, i.e. 
/3' = dp/dtp; j GJ = p GJ c is the Goldreich- Julian current 
density in the polar cap for the corotating magnetosphere. 
At the surface of the NS, where the potential drop is zero 
and plasma corotates with the NS, j GJ corresponds to the 
local GJ current density. 



2.4 Goldreich-Julian charge density in the polar 
cap for differentially rotating magnetosphere 

The Goldreich-Julian (GJ) charge density is the charge den- 
sity which supports the force-free electric field: 



Pgj 



47T 



(16) 



The GJ charge density in points along a magnetic field line 
rotating with an angular velocity different from the angular 
velocity of the NS will be different from values given by the 
eq. (|14[) . Substituting expression for the force- free magnetic 
field ([2| into eq. (fTfjjl we get 



(/3A$ + /3'(V*) 2 ). 



(17) 



We see that the GJ charge density depends not only on the 
angular velocity of the field line rotation (the first term in 
eq. HZ}), but also on the angular velocity profile (the second 
term in eq. (117)1 '). 

Near the NS the magnetic field is essentially dipolar. 
The magnetic flux function "if for dipolar magnetic field is 



(18) 



Substituting this expression into equation (|17|l we get 
M 1 „ 



Pgj 



• 2 a 

/32(3cos 2 (9- 1) (3cos 2 6» + l) ) . (19) 



In the polar cap of pulsar cos# ~ 1 and /V( r ^Lc) 3 — B/2. 
Recalling expression for the magnetic flux function for dipole 
magnetic field (|18[1 we get for the local GJ charge density in 
the polar cap of pulsar 



Pgj 



Pgj (0 + /9V). 



(20) 



3 ACCELERATING POTENTIAL 

In our approximation any current density distribution in 
force-free magnetosphere of an aligned rotator has the form 
given by eq. (|15[l . The current density depends on the an- 
gular velocity of the magnetic field lines fl(ip), which for a 
given field line depends on the total potential drop along 



that line via equation (|8]). The potential drop in the accel- 
eration zone can not exceed the vacuum potential drop, i.e. 
V is limited by inequality ((9}. 

So, if we wish to find the accelerating potential which 
supports a force-free configuration of the magnetosphere for 
a given form of the current density distributiorjf] in the po- 
lar cap we do the following. We equate the expression for 
the current density distribution to the general expression 
for the current density (|15p . then we express f3(ij}) in terms 
of V(ip) by means of equation (|8]), and obtain thus an equa- 
tion for the electric potential V which supports a force-free 
magnetosphere configuration with the desired current den- 
sity distribution. If solutions of the obtained equation fulfil 
limitation ©, such configuration is admitted, if not - such 
current density could not flow in force-free magnetosphere 
of an aligned pulsar. Currently there is no detailed model 
for non- stationary polar cap cascades from which we could 
deduce reasonable shapes for current density distribution. 
Therefore, we try to set constrains on the current density 
assuming linear dependence of the current density on ip. 

In differentially rotating magnetosphere there are 
two characteristic current densities. The first one is the 
Goldreich-Julian current density for the corotating magne- 
tosphere j GJ . ^ corresponds to the actual Goldreich-Julian 
current density in the magnetosphere at the NS surface, 
where differential rotation is not yet build up. The second 
characteristic current density is the actual Goldreich-Julian 
current density j GJ in points above the acceleration zone 
where the magnetosphere is already force-free and the final 
form of differential rotation is established; in the polar cap 
j'gj is given by formula (|20[l . For magnetosphere with strong 
differential rotation the current densities j'qj and j GJ differ 
significantly. In this section we consider both cases, namely, 
when the current density distribution is normalised to j'gj 
and when it is normalised to j GJ . 

3.1 Outflow with the current density being a 
constant fraction of the actual 
Goldreich-Julian current density 

For non-stationary cascades the physics would be deter- 
mined by the response of the cascade zone to the inflowing 
particles and MHD waves coming from the magnetosphere. 
However, the accelerating electric field depends on the de- 
viation of the charge density from the local value of the GJ 
charge density. So, the first naive guess would be that the 
preferable state of the cascade zone would be the state when 
(on average) the current density is equal to the GJ current 
density j GJ : 

jW=jo,W=j GJ G3 + /8V). (21) 

Equating this formula to the general expression for the cur- 
rent density (|15|l and substituting for /3 expression (|8]) , after 
algebraical transformation we get the equation for the ac- 
celerating electric potential in the polar cap of pulsar 

1 4- V 

V" = -2 —-j— . (22) 

W 

We set the boundary conditions for V at the edge of the 
polar cap. As the boundary conditions we can use the value 

3 guessed from a model for the polar cap cascades, for example 
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of the normalised angular velocity at the edge of the polar 
cap and the value of the electric potential there 



1 + V'(1) 
V(l) 



/3 P c 
V 



(23) 
(24) 



Solution of equation Q22J) satisfying the boundary condi- 
tions (pi) . p3l) is 



VWO = Vo + (1 - V) 1 - 



An 



(25) 



We see that unless /3 pc = the potential has singularity at 
the rotation axis, and, hence, such configuration can not be 
realised in force-free magnetosphere of a pulsar. The con- 
dition @ is violated - the potential difference exceeds the 
vacuum potential drop. 

If /3 P c ~ 0, the potential is V — Vo + 1 — ij) and from 
eq. ([SI we have /3(V0 = 0. Substituting this into eq. (fl5|) 
we get for the current density = 0. So, the case with 
Ppc = is degenerated, as there is no poloidal current in the 
magnetosphere, it corresponds to the vacuum solution. 

Let us consider now a more general form for the current 
density distribution 



(26) 



where A is a constant. In that case for the accelerating elec- 
tric potential in the polar cap of pulsar we have the equation 



V" = 2(1 + 1/') 



■ A- V 



V [V + 2(A - 1)] 



(27) 



For the same boundary conditions (|24l) . (123 solution of this 
equation is 



vW = Vo+1 ^ + Pe^±J) log 



4>{2A- 



i/) + 2(A-l) 



(28) 



This solution is valid for A ^ 1, 1/2. There is the same 
problem with the electric potential in that solution. Namely, 
unless /3 pc = the potential V is singulajf] at the rotation 
axis. The case with A = 1/2 is also degenerated, because in 
that case the solution for the electric potential is V(ip) = 
Vo + 1 — ip what yield the current density j(ip) = 0. 

We see that solutions with the current density being a 
constant fraction of the actual GJ current density are not 
allowed, except a trivial degenerated case, corresponding to 
no net particle flow. The naive physical picture does not 
work and the current density in the magnetosphere in terms 
of the actual GJ current density must vary across the polar 
cap. On the other hand, the GJ current density is itself a 
variable function across the polar cap, it changes also with 
the altitude within the acceleration zone, when the potential 
drop increases until it reaches its final value. So, we find it 
more convenient to consider the current density in terms of 
the corotational GJ current density. 



4 the singularity arises because V"(Q) goes to infinity unless 1 - 
V'(0) = /3(0) is zero, as it follows from equation ]27t 



3.2 Outflow with the current density being a 

linear function of the magnetic flux in terms 
of the corotational Goldreich-Julian current 
density 

In models with the space charge limited flow (SCLF), when 
charged particles can fre ely escape from the NS surface (e.g. 
IScharlemann et al.ll 19781 ) , the charge density at the NS sur- 
face is always equal to the local GJ charge density there, 
(P — PGj)|r=fl N s' F° r SCLF the actual current density in 
the polar cap could be less than j GJ if acceleration of the 
particles is periodically blocked in the non-stationary cas- 
cades. The current density could be greater than j'qj if there 
is an inflow of particles having opposite charge to that of the 
GJ charg e density from th e magnetosphere into the cascade 
zone fe.g. lLvubarskIilll992h . Therefore, an expression for the 
current density in terms of the corotational GJ current den- 
sity Jq j would be more informative from the point of view 
of the cascade physics. 

Let us consider the case when the current density dis- 
tribution in the polar cap of pulsar has the form 



3 =3gz (aip + b). 



(29) 



where a, b are constants. The Michel current density distri- 
bution is a particular case of this formula and corresponds 
to the values of parameters a — —l,b— 1. The equation for 
the electric potential for this current density is 



V" 



<# + &-(! + V')(l-^) 
^(2-V) 



(30) 



Solution of the equation (|30[) satisfying the boundary con- 
ditions (HP, is 



V(V0 = VO+(l + a)(l-V)+ 



+ 



i log [(2 -VO-^-^ 



(31) 



We see that the potential is non singular at the rotation axis, 
if ,0p C = a + 2b. So, the admitted solution for the electric 
potential is 

V(xp) = V0 + (1 + o)(l - V) - 2(a + b) log(2 - V) • (32) 

In the rest of the paper we will use for the electric potential 
expression (|32[1 . We will analyse physical properties of force- 
free magnetosphere configurations when the electric poten- 
tial in the acceleration zone of the polar cap has that form. 



4 PROPERTIES OF ADMITTED 
CONFIGURATIONS 

4.1 Admitted current density 

The potential drop in the polar cap of pulsar is limited by the 
vacuum potential drop. In our notations this limit is formu- 
lated as inequality (|9}. Parameters a, b from the expression 
for the electric current (|29[) enters into the formula for the 
electric potential (|32[) . Imposing limitation ((9)1 we get the 
admitted range for these parameters in the force-free mag- 
netosphere. In Appendix|A]we do such analysis and find the 
region in the plane (a, b) which is admitted by the require- 
ment (|9]). This region is shown as a grey area in Fig. IA1I 
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From Fig. IA1I it is evident that for the most of the admit- 
ted values of parameters a, b the current density has differ- 
ent signs in different parts of the polar cap. There is also 
a region where the values of parameters correspond to the 
current density distributions having the same sign as the GJ 
charge density in the whole polar cap. 

The physics of the polar cap cascades impose additional 
limitations on the current density and accelerating electric 
potential distribution in the polar cap. There is now no de- 
tailed theory of non-stationary polar cap cascades. There- 
fore, in setting constrains on the current density distribution 
we should use some simple assumptions about the possible 
current density. There is a preferable direction for the accel- 
erating electric field in the polar cap. The direction of this 
field in such that it accelerates charged particles having the 
same sign as the GJ charge density away from the star. It 
is natural to assume that the average current in the polar 
cap cascade should flow in the same direction. The average 
current could flow in the opposite direction only if the ac- 
celerating electric field is screened. In order to screen the 
accelerating field a sufficient amount particles of the same 
sign as the accelerated ones should come from the magneto- 
sphere and penetrate the accelerating potential drop. These 
particles, however, are itself produced in the polar cap cas- 
cade. They must be accelerated somewhere in the magneto- 
sphere back to the NS up to the energy comparable with the 
energy the primary particles get in the polar cap cascade. 
Even if the problem of particle acceleration back to the NS 
could be solved, screening of the electric field will interrupt 
the particle creation, and, hence, there will be not enough 
particles in the magnetosphere which could screen the elec- 
tric filed the next time. Although the real physics is more 
complicated and is not yet fully understood, the case of the 
unidirectional current in the polar cap is worth of detailed 
investigation as it is "the most natural" from the point of 
view of the polar cap cascade physics. In the following we 
will call the current being of the same sign as the GJ charge 
density as "positive" and the current being of the opposite 
sign to the GJ charge density as "negative". 

The linear current density distribution (|29p will be al- 
ways positive if 

6 ^ max(— o, 0) . (33) 

Only a subset of the admitted values of a, b corresponds to 
the positive current density distribution. Such values of the 
parameters a, b are inside the triangle-like region shown in 
Figs. [21 El HI We see that a rather wide variety of positive 
current density distributions are admitted in the force-free 
magnetosphere: current density distributions being constant 
across the polar cap of pulsar are admitted, as well as cur- 
rent densities decreasing or increasing toward the polar cap 
boundary. So, the current density in the force-free magneto- 
sphere could deviate strongly from the classical Michel cur- 
rent density, corresponding to the point a = —1,6= 1. The 
price for this freedom is the presence of a non-zero accel- 
erating electric potential in the polar cap. If the price for a 
particular current density distribution is too hight, i.e. if the 
potential drop is too large, only magnetosphere of pulsars 
close to the "death line" could admit such current density. 
Let us now consider the distribution of the potential drop 
in the parameter space (a, b). 
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Figure 2. Maximum potential drop across the polar cap. The 
dotted lines show contours of AV ma x- Contours for AVmax = 
0.05, 0.1, 0.2, 0.5, 0.8 are shown. Labels on the lines correspond to 
the values of AVmax- The line corresponding to AV max = 0.05 is 
not labelled. 

4.2 Electric potential 

We emphasised already that the shape of the function V(ip) 
is very important for the resulting current density distribu- 
tion. However, as we do not understand in detail the physics 
of non- stationary cascades, we cannot judge whether a par- 
ticular form of V(tp) is admitted by the cascade physics or 
not. On the other hand, in young pulsars the average poten- 
tial drop could not be very large, because already a small 
fraction of the vacuum potential drop would be sufficient for 
massive pair creation and screening of the accelerating elec- 
tric field. So, currently we could judge about reasonableness 
of a particular current density distribution only from the 
maximum value of the potential drop it requires. The elec- 
tric potential given by eq. (|32p is known up to the additive 
constant Vo, which is the value of the accelerating potential 
at the polar cap boundary. Vo and thus the actually poten- 
tial drop in the accelerating zone can not be inferred from 
the magnetosphere physics and is set by the physics of the 
polar cap cascades. The only thing we could say about the 
actual potential drop in the acceleration zone along field 
lines is that its absolute value is not smaller than the ab- 
solute value of the maximum potential drop of V(ip) across 
the polar cap. 

Let us now consider possible values of the maximum po- 
tential drop across the polar cap of pulsar. If the potential 
is a monotone function of ip in the polar cap, the maximum 
potential drop is the drop between the rotation axis and the 
polar cap boundary. If the potential as a function of ip has a 
minimum inside the polar cap, the maximum potential drop 
will be either between the axis and the minimum point, or 
between the edge and the minimum point. We analyse this 
issue in details in Appendix [B] In Fig. [5] the contour map 
of the maximum potential drop in the plane (a, b) is shown. 
The line given by eq. (|B1|I is the line where for fixed a (or b) 
the smallest value of the potential drop across the polar cap 
is achieved. From this plot it is evident that even if the po- 
tential drop in the polar cap is rather moderate, of the order 
of ~ 10 per cents, there are force-free magnetosphere con- 
figurations with the current density distribution deviating 
significantly from the Michel current density distribution. 
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Figure 3. Ratio of the actual current density to the Goldreich- 
Julian current density t(l) at the polar cap boundary, where the 
minimum value of this ratio is achieved, see text. The dotted lines 
show contours of t(l) 



So, even for young pulsars there may be some flexibility in 
the current density distribution admitted by the force-free 
electrodynamics. 

Note that force-free magnetosphere impose different 
constraints on pulsars in aligned /i ■ CI > and anti-aligned 
configuration fi ■ ft < configuration (p u lsar a nd antipul- 
sar in terms of iRuderman fc Sutherland! (|l975h ) . For pul- 
sars the accelerating potential is positive, i.e. it increases 
from the surface of the NS toward the force-free zone above 
the pair formation front. In case of antipulsar the poten- 
tial is negative, it decreases toward the pair formation front, 
because positive charges are accelerated. Equations for the 
current density (|15[1 , (|29[) we used to derive the equation for 
the electric potential (|30p contain the expression for the GJ 
charge density as a factor, and, hence, the resulting expres- 
sion for the electric potential is the same for each sign of 
the GJ current density. So, for pulsars there is a minimum 
in the accelerating potential distribution, for antipulsar the 
distribution of the accelerating electric potential has a maxi- 
mum. Mathematically this results from different signs of the 
integration constant Vq. 



4.3 Angular velocity 

The normalised angular velocity of the open magnetic field 
lines in the force free magnetosphere with linear current den- 
sity distribution (|29p is given by 



26 + atp 
2- V 



(34) 



For admitted current densities it grows with increasing of ip, 
because the first derivative df3/dip for the admitted values 
of a, b is always non-negative. So, the angular velocity either 
increase toward the polar cap boundary or remains constant 
over the cap if a = —6. The latter case includes also the 
Michel solution. The minimum value of j3 



(35) 



is achieved at the rotation axis, where if) = 0, and the max- 
imum value 



Figure 4. Spindown rate in terms of the Michel spindown rate. 
Dotted lines show contours of w. Label on the lines correspond 
to the values of w. 



at the polar cap boundary, where ifr = 1. So, the open field 
lines can rotate slower, as well as faster that the NS, but the 
lines near the polar cap boundary rotate not slower than the 
lines near the rotation axis. 



4.4 Goldreich-Julian current density 

Expression for the GJ current density in the polar cap can 
be obtained by substitution of the expression (|34[) for /3 into 
equation (|20|) for the GJ current density. We get 



.o 4b + aip(4-ip) 



(V- - 2) 2 



(37) 



For the admitted values of the parameters a, b the derivative 
djGj/dtp is always non-negative and, hence, the GJ current 
density either increases toward the polar cap boundary, or 
remains constant when a — —b. The actual current density, 
however, could decrease as well as increase toward the polar 
cap edge. 

For the charge separated flow the deviation of the cur- 
rent density from the GJ current density generate an ac- 
celerating or a decelerating electric field when j < j G j or 
j > Jg.t correspondingly. Although in non-stationary cas- 
cades the particle flow would be not charge separated, the 
ratio of the actual current density to the GJ current den- 
sity may give some clues about cascade states required by a 
particular magnetosphere configuration. This ratio is given 
by 



_ W>-2) 2 (6 + ay>) 



j G j(i>) tM/>(4-V) + 46 



(38) 



For each admitted configuration the current density is equal 
to the GJ current density at the rotation axis. For the ad- 
mitted values of the parameters a, b the derivative di/dtj} 
is always positive, and, hence, the current density in terms 
of the GJ current density decreases toward the polar cap 
boundary. So, except the rotation axis the current density 
in the polar cap is always less than the GJ current density. 
The relative deviation of the actual current density from the 
GJ current density is maximal at the polar cap boundary 



3 a 



2b + a, 



(36) 



t(l) = 



a + b 
3a + 46 



(39) 
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Figure 5. Electric potential in the polar cap of pulsar as a func- 
tion of the normalised flux function xp for magnctosphere config- 
urations with a constant current density across the cap. In all 
cases Vq is set to zero. Numbers near the lines correspond to the 
following values of b: 1 — b = 0; 2 — b = .5; 3 — b = f> ma x/2; 4 
— 6 = 1; 5 — b = femax- The line corresponding to the minimum 
potential drop across the cap is shown by the thick solid line (the 
line 3). 



Its maximum value t max = 1/3 this ratio achieves when 
6 = 0. Its minimum value t m i n = it achieves when a = —b, 
what includes also the case of the Michel's current density 
distribution. The contours of t(l) are shown in Fig. [3] 

4.5 Spin-down rate and the total energy of 

electromagnetic field in the magnetosphere 

In our notations the spindown rate of an aligned rotator is 
(cf. eq. (60) in Paper I) 



W = W, 



/•*pc 

Jd / 
Jo 



(40) 



where Wmd is the magnetodipolar energy losses defined as 



W md = 



4c 3 



(41) 



Substituting expression for the poloidal current (|12[) and 
using the normalised flux function ip we get 



/3 2 (ip)ip(2 -ip)dtp. 



(42) 



Expression for the spindown rate in the Michel solution 

(43) 



W M = !*p C W md 



differs from the spindown rate obtained in the numerical 
simulations of the corotating aligned rotator magnetosphere 
with by a constant factor. However, it has very similar de- 
pendence on the size of the corotating zone xo (cf. equa- 
tions (62), (63) in Paper I). As our solutions are obtained in 
split-monopole approximation, they should differ from the 
real solution approximately in the same way as the Michel 
solution does. Because of this it would be more appropriate 
to normalise the spindown rate to the spindown rate in the 
Michel split-monopole solution. Doing so we will be able to 
study the effect of differential rotation on the energy losses 
separately from the dependence of the spindown rate on the 
size of the corotating zone. 



Figure 6. Normalised angular velocity of the open magnetic field 
lines as a function of the normalised flux function ip for magneto- 
sphere configurations with a constant current density across the 
cap. Labelling of the curves is the same as in Fig. \E\ 



For the normalised spindown rate in the considered case 
of linear current density we get 



W 

W M 



: 4a 2 (3 log 2-2)+ 

+ 3a6(8 log 2 - 5) + 6b 2 (2 log 2 - 1) . (44) 

In Fig. [4] contour lines of w are shown in the domain of 
admitted values for parameters a, b. We see that spindown 
rate can vary significantly, from zero to the value exceeding 
the Michel's energy losses by a factor of « 6. It increases with 
increasing of the values of the parameters a, b and decreases 
with decreasing of that values. It is due to increasing or 
decreasing of the total poloidal current in the magnetosphere 
correspondingly. 

The total energy of the magnetosphere could be es- 
timated from the split-monopole solution. Using the for- 
mula (|C7P derived in Appendix [C] we have for the total 
energy of the electromagnetic field 



W~ W P o1 + (R-Rns)W, 



(45) 



where W po i is the total energy of the poloidal magnetic field 
and R is the radius of the magnetosphere. The first term in 
our approximation is the same for all magnetosphere config- 
urations, the difference in the total energy arises from the 
second term. Hence, the contours of the constant total en- 
ergy in the plane (a, b) have the same form as the contours 
of the spindown rate W shown in Fig. [4] So, the total energy 
of the magnetosphere increases with increasing of parame- 
ters a, b, i.e. it increases with the increase of the poloidal 
current . 



4.6 Example solutions 

As examples we consider here the properties of two partic- 
ular solutions in details. We chose these solution because 
either their current density or the potential drop seem to 
correspond to "natural" states of the polar cap cascades. 
Although we do not claim that one of the solutions is re- 
alised as a real pulsar configuration, but knowledge of their 
properties may be helpful in understanding of the physical 
conditions the polar cap cascades should adjust to. 
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Figure 7. Current density as a function of the normalised flux 
function ip for magnetosphere configurations with a constant cur- 
rent density across the cap. Labelling of the curves is the same 
as in Fig. \E\ By the thick grey line the ratio of the actual current 
density to the GJ current density i(i/>) is shown. For this case it 
is the same for all solutions. 

4-6.1 Configurations with constant current density 

At first we consider the case when the current density is con- 
stant over the polar cap, i.e. a — and j = b j'gj- Constant 
density distribution would be produced by cascades in their 
"natural" state, if the current adjustment proceeds locally, 
without strong influence from the current along adjacent 
field lines. The electric potential in that case is 

V c {ip) = Vo-\-l-ij)-2b\og{2-i))). (46) 

This potential has the following properties (see Fig. [5] where 
V(ip) is shown for several values of b assuming for the sake 
of simplicity Vo = 0): 

• the admitted values of the current density in the po- 
lar cap of pulsar are within interval [0, b max ], where b max = 
1/ log 2 ~ 1.443. 

• if < b < b max /2 ~ 0.721 the value of the electric 
potential at the rotation axis V c (0) is larger than that value 
at the polar cap edge V C {1), V c {0) > V C {1) 

• if b max /2 < b < b max the value of the electric potential 
at the rotation axis V c (0) is smaller than that value at the 
polar cap edge V c (l), V c (0) < V c (l). 

• if0<b<l/2orl<b< b max the potential is a 
monotone function of ip; if 1/2 < b < b m&K it has a minimum. 

• in the point b = b max /2 the maximum potential drop 
across the polar cap reaches its minimum value AKn ax = 
0.086. 

The reason for such behaviour of the potential is easy to un- 
derstand from the Fig. lBll in Appendix[B] The critical points 
where V(ip) changes its behaviour are the points where the 
line a = intersects the boundaries of the regions I,II,II, 
and IV. 

The angular velocity of the open magnetic field lines is 

= ^ (47) 

Distribution of the corresponding angular velocity is shown 
in Fig. \§\ For b > 1 the angular velocity of rotation of all 
open magnetic field lines is larger than the angular velocity 
of the NS. For b < 1/2 all magnetic field lines rotate slower 
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Figure 8. Spindown rate of an aligned rotator normalised to 
the spindown rate in the Michel solution for magnetosphere con- 
figurations with a constant current density across the cap as a 
function of the current density in the polar cap (parameter b). 

that the NS. For 1/2 < b < 1 some open filed lines near the 
rotation axis rotates slower that the NS, the other open field 
lines rotates faster than the NS. 

The current density distribution in terms of the GJ cur- 
rent density is 

t c (V0 = i(2-^) 2 , (48) 

and it doesn't depend on the value of the parameter b. The 
current density is always sub-Goldreich- Julian, except the 
rotation axis, where it is equal to the GJ current density. 

The normalised spindown rate for the considered case 
has simple quadratic dependence on the current density 

w c = 6(log4- l)b 2 . (49) 

This dependence is shown in Fig. [8] The energy losses in 
configuration with a constant current density could not be 
higher than ~ 4.82 of the energy losses in the corresponding 
Michel solution. 

The case b = 1 is worth to mention separately, as it is 
"the most natural" state for the space charge limited par- 
ticle flow, for which the current density at the surface on 
the NS is equal to the corotational GJ current density. In 
Fig s - El El 0the lines corresponding to this case are labelled 
with "3". The maximum potential drop for the configura- 
tion with the current density distribution being equal to the 
corotational GJ current density is AV max = 0.386 and the 
angular velocity of the open field lines varies from 1 at the 
rotation axis to 2 at the polar cap boundary. 

4-6.2 Configurations with the smallest potential drops 

As the next example we consider the case when the maxi- 
mum potential drop across the polar cap for a fixed value of 
either a or b is minimal. The points corresponding to such 
values of parameters are show in Figs. [21 El HI by the thick 
grey line. Equation for this line in the plane (a, 6) is de- 
rived in Appendix [Bj equation (|B1[1 . In some sense this is 
an optimal configuration for the cascade zone, because for 
a fixed value of the current density at some magnetic field 
line such configuration requires the smallest potential drop 
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Figure 9. Electric potential in the polar cap of pulsar as a func- 
tion of the normalised flux function rjj for magnetosphere con- 
figurations with the smallest potential drop across the cap. In all 
cases Vb is set to zero. Numbered lines correspond to the following 
values of a: 1 — a = — 1 (Michel's solution); 2 — a = (solution 
with a constant current density); 3 — a = 1; 4 — a = 2;5 — 
a = l/(Zog4 - 1). 




Figure 10. Normalised angular velocity of the open magnetic 
field lines as a function of the normalised flux function ip for 
magnetosphere configurations with the smallest potential drops 
across the cap. Labelling of the curves is the same as in Fig. [9] 



among the other admitted configurations. The accelerating 
potential for the considered class of configurations is 



V s (iP) = V -(a+l)t4> + lo. 



1 - -r 



(50) 



The potential is shown as a function of tjj in Fig. [9] for sev- 
eral particular cases assuming for the sake of simplicity zero 
potential drop at the polar cap boundary. The potential has 
always a minimum at the point 



log 2 



0.557, 



(51) 



the position of this minimum does not depend of the values 
of a, 6. The minimum value of the maximal potential drop 
across the cap, min(A Vmax) = 0, is achieved at the left end of 
the grey line, at the point (a = — 1, 6=1) corresponding to 
the Michel solution. The maximum potential drop across the 
gap for this class of configurations, max(AVmax) = 0.309, 



Figure 11. Current density as a function of the normalised flux 
function tjj for magnetosphere configurations with the smallest 
potential drops across the cap. Labelling of the curves is the same 
as in Fig. [9] 








Figure 12. Ratio of the actual current density to the GJ cur- 
rent density i as a function of the normalised flux function ip for 
magnetosphere configurations with the smallest potential drops 
across the cap. Labelling of the curves is the same as in Fig|9] 



is achieved at the right end of the grey line, at the point 
(a = l/(log4-l), 6 = 0). 

The angular velocity of the open field lines is 



a + l 



(2-^) 2 



(52) 



Distribution of /3 s (tp) is shown in Fig. 1101 Before the min- 
imum point f/'min P is not greater than 1, after that point 
(3 is not smaller than 1. With increasing of the maximum 
potential drop the variation of the angular velocity across 
the polar cap becomes larger. 

The current density distribution in the considered case 
has the form 



f(ip) =ait> + 



q(l-log4) + l 
log 4 



(53) 



That distributions is shown in Fig. 1111 All curves pass 
through the point ip = 1 — 1/ log 4, where the current density 

kj"(i?)=i&j/i°g4. 

The current density distribution in terms of the 
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Figure 13. Spindown rate of an aligned rotator normalised to 
the spindown rate in the Michel solution for magnetosphere con- 
figurations with the smallest potential drops across the polar cap 
as a function of the parameter a. 

Goldreich- Julian current density is 

,s M = (^-2) 2 [(^-l) a log4 + a + l] 
y ' i/>(4 - V)alog4 + 4a(l - log 4) +4 ' V ; 

It monotonically decreases from 1 at the rotation axis to its 
minimum value at the PC boundary. This minimum value is 
in the range [0, 1/3], the lowest value corresponds to the left 
end of the grey curve (the Michel solution) , the upper value 
corresponds to the right end of the grey line. The outflow is 
sub-GJ everywhere except the rotation axis. 

The normalised spindown rate for the considered case 

is 

w s = 3_ { \2 log 2 2 - 3(1 - log 2)1 a- 

2 log 2 2 U y ;J 

-3(2 -3 log 2)a- 3(1 -2 log 2)} . (55) 

It is shown as a function of the parameter a in Fig. 1131 
We see that for the considered configuration the spindown 
rate as a function of the parameter a increases more slowly 
than the spindown rate for configurations with a constant 
current density as a function of b, cf. Fig[S] The energy losses 
could not be higher than ~ 2.13 of the energy losses in the 
corresponding Michel solution. 



5 DISCUSSION 

The main aim of this paper was to study the range of admit- 
ted current density distributions in force-free magnetosphere 
of an aligned rotator. Taking into account that this subject 
was not studied in details before, the linear model used in 
this work is in our opinion an adequate approach to the 
problem. The knowledge of the magnetosphere behaviour in 
response to different potential drops in the polar cap would 
be very useful for future modelling of non-stationary polar 
cap cascades. This formalism could be used as a tool allowing 
to judge quickly whether a particular model of the polar cap 
cascades is compatible with the force-free magnetosphere or 
not. It may also give a clue how the magnetosphere would re- 
spond to a particular current density distributions obtained 
at some step in course of numerical solution. Although the 



analytical model presented here needs to be refined in nu- 
merical simulations, the presence of some analytical rela- 
tions should be very handy for numerical cascade modelling. 

We have considered here a simple case when the current 
density in the polar cap of pulsar is a linear function of the 
magnetic flux. However, the generalisation of the model to 
the case of a more complicated shape of the current density 
distribution is straightforward. One should proceed with the 
steps described at the beginning of the Section [3] for the de- 
sired form of the current density distribution. The resulting 
equation for the electric potential will be an ordinary differ- 
ential equation, and numerical solution of such equation for 
any given current density will not be a problem. 

The main conclusion we would like to draw from the 
presented results is that even for rather moderate potential 
drop in the acceleration zone the current density distribu- 
tion can deviate significantly from the "canonical" Michel 
distribution, which is basically preserved for dipole geome- 
try if all field lines corotates with the NS (see Paper I). In 
particular, a magnetosphere configuration with the constant 
current density at the level of 73 per cents of the Goldreich- 
Julian current density at the NS surface would require a 
potential drop in the acceleration zone of the order of 10 
per cents of the vacuum potential drop. For time-dependent 
cascades this may be realised even for young pulsars. We 
should note however that for young pulsars the potential 
drop of the order of 10 per cents of the vacuum drop could 
cause overheating of the polar cap by the particles acceler - 
ated toward the NS (e.g. iHarding fc Muslimovll200ll . 12051 ). 
In that sense such potential drop may be too large for young 
pulsars. On the other hand, without knowledge of the dy- 
namics of non-stationary cascades it is in our opinion too 
preliminary to exclude the possibility of such configurations 
for young pulsars, as in non-stationary regime the heating 
of t he cap may be not s o strong as in stationary cascades 
f see iLevinson et afll2005t) . 

We used split-monopole approximation for the poloidal 
current density distribution in the magnetosphere, which 
would produce an accurate results only for configurations 
with very small size of the corotating zone, less than the 
size of the light cylinder xo <C l/(3(ty pc ), see Fig. [T] For the 
(most interesting) case when these sizes are comparable, re- 
sults obtained in this work could be considered only as a zero 
approximation to the real problem. We should point out an 
important modification introduced by the dipole geometry 
of the magnetic field. For the dipole geometry there will be 
some magnetic field lines which bent to the equatorial plane 
at the light cylinder. For these field lines the second terrrQ 
on the r.h.s. of equation Q will be negative, and in order to 
get positive current density along these lines a more steep 
dependence of (5 on i\> would be necessary. As a result the po- 
tential drop in the dipole geometry would be higher than it 
is obtained in our approximation. Figuratively speaking, in 
our model we could correct only for the decrease of the elec- 
tric current density toward the polar cap boundary present 
in the Michel solution. On the other hand, if /3($ pc ) > L 
the size of the corotating zone can be less as well as greatejf] 



5 which in the cylindrical coordinates (w, (j>, z) is Ifid-^^ 

6 for the closed magnetic field lines the angular velocity is fl and 
they are still inside their light cylinder; the adjustment of the 
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then the light cylinder radius at the last open field line. In 
the latter case there should be less magnetic field lines which 
bent toward the equatorial plane, than in the first case, cf. 
Fig. [T] cases I and II. Hence, the correction introduced by 
the dipolar field geometry for some subset of our solutions 
would be non-monotonic as the size of the corotating zone xo 
increases. So, there is still possibility that a moderate poten- 
tial drop could allow a large variety of the current densities, 
although this issue needs careful investigation. 

In this paper we ignored the electrodynamics of the po- 
lar cap zone. Although without a theory of time-depending 
cascades we cannot put more limits on the electric potential 
than the one we used in Section 12.21 there is an additional 
limitation coming from the basic electrodynamics. Namely, 
the accelerating potential near conducting walls, which the 
current sheet at the polar cap edge is believed to be, should 
approach zero. However, we could speculate that there is a 
thin non-force-free zone at the edge of the polar cap where 
the adjustment of the potential happens. In other words, the 
return current and the actually nearly equipotential region 
may occupy not the whole non-force-free zone. Because of 
this that limitation would not eventually restrict our solu- 
tions. 

At the end we would like to discuss briefly the issue with 
the pulsar breaking index. If the inner pulsar magnetosphere 
is force-free, then the spindown rate of an aligned pulsar as 
a function of angular velocity will deviate from the power 
law W oc Q 4 if the size of the corotating zone or/and the 
distribution f3(ip) change with time. The assumption that 
these "parameters" are time dependent seems to us to be 
natural, because with the ageing of pulsar the conditions in 
the polar cap cascade zone change and the magnetosphere 
should adjust to these new conditions. In the frame of our 
model we could make some simple estimations how the brak- 
ing index of pulsar is affected by the changes of these two 
"parameters" . 

As an example we consider the case when the pulsar 
magnetosphere evolves through a set of configurations with 
a constant current density. The spindown rate for such con- 
figurations is 

W tx STT^pcb 2 ~ O. 4 xo V , (56) 

here we estimate type assuming dipole field in the corotation 
zone. If b and/or xo are functions of time, the spindown rate 
will be different from the spindown of a dipole in vacuum. If 
at some moment the dependence of the size of the corotating 
zone and the current density on Q could be approximated 
as xo oc and 60 oc fi^ correspondingly, the braking index 
of pulsar measured at that time is 

n = 3 - 2£ + 2C . (57) 

We see, that the deviation of the braking index from the 
"canonical" value being equal to 3 would be by a factor 
of two larger than the dependence of b and xo on fi. The 
braking index could be smaller as well as larger than 3, de- 
pending on the sign of the expression ( — £. For instance, if in 
an old pulsar the potential drop increases and, as the con- 
sequence, the current density decreases, £ is positive, and 

angular velocity occurs in the current sheet, which is for sure a 
non-force-free domain. 



the braking index could be greater than 3. We note, that 
there are ev idences for such values of the braking index for 
old pulsars (lArzoumanian et al.ll2002T ). On the other hand, 
if £0 decreases with the time, as it was proposed in Paper I 
for young pulsars, the braking index would be less than 3. 
However even in this simple picture different trends may 
be possible. In reality the evolution of the pulsar magne- 
tosphere will be more complicated and more steep as well 
as more gradual dependence of the braking index on the 
current density would be possible. This would results in a 
rather wide range of possible values of pulsar braking index. 
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Figure Al. Value of parameters a and b admitted by the re- 
quirement that the maximum potential drop across the polar cap 
is less than the vacuum potential drop. Dotted lines 1 and 1' 
show the boundary of the region where the potential has an ex- 
tremum in the polar cap. Line 2 corresponds to AVio = 1, and 
2' - to AVio = — !• Line 3 show points when |AVi e | = 1, and 
4 - |AVbe| = 1- Line 2' corresponds to b = —a. The resulting 
admitted region is shown by the grey colour. See text for the 
explanations. 
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APPENDIX A: 
DENSITY 



ADMITTED CURRENT 



The second derivative of V(tf>) in respect to ip is 



a + b 
(2-V) 2 



(Al) 



For fixed a and 6 the second derivative never changes the 
sign, hence, V(tp) has a single extremum. If a + b > 0, V(ip) 
has a single minimum; in the plane (a, b) shown in Fig. IA1I 
these points lie left to the line 5. If a + b < 0, V(ip) has a 
single maximum, such points lies to the right of the line 5 
in Fig. IA1I The potential reaches its extremum value at the 
point 



a + 1 



This point lies in the interval [0, 1] if 



l>b> 



l<b< 



2 ' 
1 - a 



for a > — 1 : 



for a < — 1 . 



(A2) 



(A3) 



In Fig. IA1I points where these conditions are satisfied lie 
between the lines 1 and 1', in the region where the angle 
between the lines is acute. The line 1 corresponds to the 
values of a, b when the extremum of the potential is reached 
at the polar cap boundary. For a, b lying on the line 1' the 
potential reaches the extremum at the rotation axis. 

If t/> x is outside of the interval [0, 1], V(tp) is a monotone 
function of %p in the polar cap of pulsar and the maximum 
potential drop is the potential drop between the edge of the 



polar cap and the rotation axis. In that case the condition ((9]) 
takes the form 



I AVio | «S 1. 



(A4) 



where AVio = V(l) — V(0) is the potential drop between 
the the boundary of the polar cap and the rotation axis. In 
terms of the parameters a, b it can be written as 



log 4 



«S b € 



+ a 



1 - log 4 



(A5) 



log 4 log 2 log 4 

In Fig. IA1I points satisfying this condition lie between the 
lines 2 and 2'. In the region, where the angle between the 
lines 1 and 1' is abuse (here the point Vex is outside of 
the interval [0, 1]) the lines 2 and 2' sets the boundaries for 
admitted values of the parameters a and b. 

If the electric potential reaches its extremum value in- 
side the polar cap, the maximum potential drop is achieved 
either between the extremum point and the edge of the polar 
cap, or between the extremum point and the rotation axis. 
In that case the condition ([9]) takes the form 

max(|AVi |,|AVoc|, |AVi e |) < 1 , (A6) 

where AVbe = V(0) — V{ip BX ) is the potential drop between 
the rotation axis and the point ip ex _, and AVi e = V(l) — 
V(ip ex ) is the potential drop between the polar cap boundary 
and the point ip cx . Expression for the extremum value of V 
is non-linear in respect to a, b: 

2(a + b)~ 
1 + a 

in terms of a, 6 should be evaluated nu- 



V(V>ex) = V - 1 + a + 26 - 2(o + b) log 



(A7) 



and condition (|A3 
merically. For a fixed a the derivatives of AVo c and AVi 
respect to b are 



db 

rfAV lc 
db 



= 2 log 
= 2 log 



a + b 
1 + a 
a + b 



log 4. 



(A8) 
(A9) 



If ipex G [0, 1] and condition (|A3|) is fulfilled, then d AVi e /db 
is positive and d AVo a /db is negative. So, for a fixed a when 
ipex £ [0, 1] AVoo decreases and AVi c increases with increas- 
ing of b. 

In Fig. lAll the line 3 represents points where |AVi c | = 1, 
and the line 4 - the points where |AVo e | = 1. To the right 
of the line 5 V(ip) has a minimum, and the lines 3 and 
4 represent points where AVi c = 1 and AVo c = 1 corre- 
spondingly. To the left of the line 5 V(ip) has a maximum, 
and here the lines 3 and 4 correspond to AVi c = — 1 and 
AVbe = — 1. The line 4 lies always outside of the region be- 
tween lines 1 and 1', and, hence, the absolute value of the 
potential drop between the extremum point and the rota- 
tion axis | AVbe | never achieves the vacuum potential drop 
when the extremum point is inside the interval (0, 1). 

AVie increases with increasing of b. So, below the curve 
3 to the right of the line 5, and above the curve 3 to the left 
of the line 5 |AVbo| is less than 1. Hence, when ip ex £ [0, 1] 
(a, b are in the area between the lines 1 and l'), the admitted 
values of a and b axe limited by the lines 3 and 2' to the left 
from the line 5, and by the lines 2 and 3 to the right from 
the line 5. 

Combining all discussed restrictions we get the region 
of the admitted values of parameters a and b, which is shown 
in Fig. lAll bv the grey area. 
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Figure Bl. Region of admitted values of the parameters a and 
b corresponding to current densities being of the same sign as the 
Goldreich-Julian current density. The dotted line 1 shows points 
where AVic = 0, the line 1' - points where AVbc = 0. The 
thick grey shows points where |AVbo| = |AVi e |. In the region I - 
Ay max = AVie, in the region II - AV max = AVbc in the regions 
II and IV - AV max = |AVio|. See explanations in the text. 



APPENDIX B: THE MAXIMUM POTENTIAL 
DROP ACROSS THE POLAR CAP 

In Fig. IBll the dotted lines 1 and 1' limit the regioifl in the 

parameter space (a, b) where the potential in the polar cap 
is a non-monotone function of ip and it has a minimum at 
some point i/) cx £ (0, 1), see Appendix lAl In the regions III 
and IV the potential V(ip) is a monotone function of ip. 

The potential drop AVbe = V(0) — V(tp ex ) between the 
rotation axis and the point i/> ex where V(tf>) achieves its min- 
imum value, for a fixed a decreases with increasing of b. For 
a, b at the line 1' ip e ^ is at the rotation axis, and AVbc = 
there. The potential drop AVie = V(l) — V(tp ex ) between 
the polar cap boundary and the point ip m for a fixed a in- 
creases with increasing of b. For a, b on the line 1 i/> ex is at 
the polar cap boundary and AVi c = 0. So, AVbc increases 
in the direction from the line l' to the line 1, and AVi e 
increases in the direction from the line 1 to the line l'. 

Along some line between lines V and 1 the potential 
drops AVbc and AVi e become equal. This means also that 
the potential drop AVio = V(l) — V(0) between the polar 
cap edge and the rotation axis is zero there. Equation for 
this line is obtained easily from the requirement AVio = 0: 



1 1 
+ a- 



log4 



log 4 



log 4 



(Bl) 



This line is shown in Fig. IB1I by the thick grey line. Above 
the grey line in Fig. IB1I AVi r > AVbc, and below it AVi a < 
AVbe- Hence, the line given by eq. (|B1|) is the line where the 
maximum potential drop across the polar cap achieves its 
minimum value for fixed a or b. 

Taking all this into account we conclude, that the max- 
imum potential drop across the polar cap AVmax is equal to 
the following potential drops: in the region I - to AVi e , in 
the region II - to AVbe, in regions III and IV - to |AVo c |- 



APPENDIX C: ENERGY OF THE 
ELECTROMAGNETIC FIELD IN 
SPLIT-MONOPOLE CONFIGURATION 

The energy density of the electromagnetic field in the mag- 
netosphere is 



8tt 



(E 2 



?pol 



(CI) 



For the split monopole solution when ^ = ^/ pc (l — cos 6*) the 
non-zero components of the electric and magnetic fields are 

, sin 8 



Eg 

13 r 

B 



Rlc r 2 

H sin t 



(C2) 
(C3) 
(C4) 



So, the electric field is equal to the toroidal magnetic field, 
Eg — Btf,. The total energy of the magnetosphere is then 



W = i / (Bl 

87T 



2Ej) dV 



On the other hand, the energy losses are 
[E x B] 



U = f c- 

J4tt 



47T 



dtt = 



(C5) 



(C6) 



where SI is a solid angle. Using equations (|C6|I . ()C3|I we can 
rewrite the expression (|C5|I for the total energy of electro- 
magnetic filed as 



W : 



Wu 



LC / 1 

V-Rns 



+ W(R-Rns) 



(C7) 



where R is the size of the magnetosphere (where it dynamics 
is determined by the central object). The first term repre- 
sents the total energy of the poloidal magnetic field and it is 
the same for all configurations with different potential drops. 
The second term, being the sum of the energies of the elec- 
tric field and the toroidal component of the magnetic field is 
different for different values of the accelerating potential. It 
is directly proportional to the energy losses in a particular 
configuration. 



7 regions I and II together 



